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Abstract 

We test the stability of the mean-field solution in the Nambu-Jona-Lasinio model. For stable 
solutions with respect to both the a and vr directions, we investigate effects of the mesonic loop 
corrections of 1/Nc, which correspond to the next-to-leading order in the 1/Nc expansion, on the 
high density chiral phase transition. The corrections weaken the first order phase transition and 
shift the critical chemical potential to a lower value. At Nc = 3, however, instability of the mean 
field effective potential prevents us from determining the minimum of the corrected one. 

PACS numbers: ll.30.Rd, 12.40.-y 



*matsuz a@fukuoka-edu. ac.jp 



1 



At high temperature and density quantum chromodynamics (QCD) undergoes quahtative 
changes of great physical interest. Although many works are done, aspects of the transition 
region are not under full theoretical control. It is mandatory to deepen its understanding 
in order to understand the structure of compact stars and the history of the early uni- 
verse, as well as results of ultra relativistic heavy ion experiments. With the aid of the 
progress in computer power, lattice simulations have become feasible for thermal system 
with zero or small density. In general, chiral restoration and deconfinement are expected 
to occur simultaneously One of the most important recent findings is strong correla- 
tions in the deconfined quark gluon plasma just above the critical temperature; on the one 
hand it appears as the near perfect fluidity 0, ^, ^ and on the other hand as the mesonic 
correlations 

As an approach complementary to the first-principle lattice QCD simulation, we can con- 
sider effective models. In particular, they are even indispensable at high density where lattice 
QCD is not applicable due to the signproblem. One of them is the Nambu-Jona-Lasinio 
(NJL) model. Since it was proposed [71, Isj, although this is a model of chiral symmetry 
that does not possess a confinement mechanism, this model has been widely used jol, in 
the mean field approximation, for example, for analyses of the critical 
transition on the temperature (T) - chemical potential (/i) plane ll|, 
few studies that go beyond the mean field approximation were reported. Nikolov et al. 
extended the mean field approximation to next-to- leading order in l/N^ which includes the 
meson-loop contributions. The extension preserves the symmetry property of the theory. 
They studied the meson loop effect in the case of T = yU = and concluded its importance 
due to the light pion mass. Hiifner et al. [l^ formulated thermodynamics of the NJL model 
to order 1/Nc and Zhuang et al. [l^] presented numerical results by using the formulation. 
Their analyses are mainly made for thermal system in the limit of zero /i. Although only a 
result is reported for the case that both T and are finite, one can expect from the result 
that the meson loop corrections are important also for finite /i. 

In this Letter our attention is focused on the finite-density chiral phase transition in the 
limit of zero T. We first investigate the stability of the mean field solution and for the stable 
mean field solution we further evaluate the magnitude of corrections of meson loops, that 
is, of the next-to-leading order in the l/N^ expansion, using the auxiliary field method of 
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Kashiwa and Sakaguchi which is essentially equivalent to the formulation of Hiifner et 
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al. [161. 

The Lagrangian density of the NJL model is 

G 

C = q{x){i^ - mo)q{x) + ^((g(a;)g(x))^ + {q{x)i-f5raq{x)f) , (1) 

where q{x) stands for two flavor quarks, and Nc for the number of color degrees of freedom, 
mo is the current quark mass and is the Pauli matrix. G is the coupling constant of a 
four- Fermi interactions scaled by the color number in order to use the l/N^ expansion 
later. The partition function for the NJL model reads 



y VqVqexp 









G N 

{{q{x)q{x)f + {q{x)^^,T^q{x)f) - ^{Jl{x) + J^{x)) 



(2) 



Here we introduced external fields and J^^. The last term in Eq. ([2]) is introduced for 
later convenience. ^ stands for the chemical potential of quarks. 

The partition function ([2]) has the four-Fermi interactions, accordingly we cannot carry 
out the Gaussian integration with respect to the quarks. Therefore, we introduce the auxil- 
iary fields a{x) and 7f(x) by using of the identity 

^ j d^x (a{x) + ^q{x)q{x)\ 



1 



Vcxexp 



(3) 



and the corresponding one for if. After integrating out quarks and shifting the auxiliary 
fields (cr(x) — Ja{x)) 1-^ a{x) and {ti{x) — Jt,{x)) ^ 7r(a;), we obtain the partition function 
of the auxiliary fields. 



Z = j PcrPTrexp [-A^c/[o-,vr]] , 



(4) 



with 



/[a. 



d'^x ( ^ ((7^(a;) + 7r2(x)) + ^ {a{x)J,{x) + ti{x) ■ X( 



x\ 



- Indet ((9^7^ + mo - yU74 + o-(a;) + ^Ts^"- vr(x)) . (5) 
The generating function for the connected Green's function W is defined as 



Z[J^, J^] = exp -NcW[J^, 



(6) 
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The fields Lp^j and (p-,, are also defined as 

G -5.Ux)~ G ~ n}^^^'^^'"'' G ^> 

for Lfa- and the corresponding one for (^^r- The fields represent the vacuum expectation values 
of the auxiliary fields in the external fields Ja{x) and Jt^{x). The effective action is defined 
by the Legendre transformation 

r [V?^, (pj^ = W[Ja, ^] ~ ^ y d-^^ (^^a{x)Ja{x) + 0nix) ■ Xix)^ . (8) 

Setting J(j(x) and J.j^{x) to constants leads to the effective potential 

r [cp^, fsvvi^^, 0^) , (9) 

where note that (Pa{x) and 0Tr{x) become constants v^o- and 0^,, respectively. 

In order to carry out the path integral in Eq. (j4]), we expand Eq. ([5]) around the classical 
solution (ctq, ttq): 

1 f ,4^ >4„ ' - _^ ^ I laa Ian \ I (J - ao 



h = j^d^x (^K' + ^o) + ^(^0^. + Tfo ■ X)^ (11) 
- In det(7^9^ + - /i74 + cxq + i-^^f ■ vfo) , 

4? = ^'^^'^(^ - y)^ai + tr[ra5(x, y : ao, 7ro)r^5(i/, x : ao, vfo)] , (12) 

where /a = 1 for a = a and 75X0 for a = tTq (a = 1,2,3). Here the classical solutions are 
governed by 



51 



5o{x) 
51 



5tt{x) 



= ^io-o{x) + Jaix)) - tYS{x,x : (Jo(x),7fo(x)) = , (13) 

0'=0'() ,7T=7VQ 

1 

= -^iTCo{x) + Jt,{x)) - tr [i75fS'(x, x : cro(x), 7fo(x))] = , (14) 



ilfid^ + mo - /i74 + cxo(x) + 2757" ■ 7ro(x))S'(x, ?/ : cro(x), 7ro(x)) = 5^^\x - y) , (15) 



4 



where S{x,y : ctq, 7?o) denotes the quark propagator in the external fields Ja{x). Making 
change of variables such that {(t{x) — cro{x)) t-^ a{x)/y/Wc and (vf(x) — 7ro(a;)) t-^ ^l{x)/^/Nl 
and carrying out the Gaussian integral with respect to a{x) and vf(x), we derive W as 



W^[^^^X] = /o + ^Trlnj(^) + o(^) 



(16) 



where Tr is the trace for mesonic degrees of freedom and the spacetime coordinate. Equation 
(fT6|) shows that the expansion (fTOj) introduced above turns out to be the l/N^ expansion. 
Inserting f|T6l) into ([7]), we get the relation between the classical fields and the corresponding 
fields ipcr{x) and 0tt{x): 



<^„{x) = ao{x) + 



0^{x) = TToix) + 



2N,6J^{x] 
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■(Trln/(2)) = ao{x) + 



2N,SJ^{x) 

The Legendre transformation of Eq. (fT6!l leads to 



(Trln/(2)) = 7fo(x) + 



(17) 
(18) 

(19) 



Vo(^a, ^tt) = ^'^'^^ + ~ ^ det(7^9^ + mo - /i74 + + «75^ ■ , (20) 



= ^^^Trlnj(2)(^^,^^) + o(^^^ . 



(21) 



The leading term Vq corresponds to the contribution of the mean field approximation and 
the the next-to-leading contribution Vi to the contribution of meson loops of 1/Nc order. 

Actual calculation of the effective potential 0191) is done in the momentum representation. 
The mean field part Vo is represented as 



2N 



(2^ 



Ep+\ ln(l + e-^(^''+^)) + i ln(l + e'^^^^-'^)) 



(22) 



where Ep = (p^ + -|- 'f^Y^'^ and M = mo -|- (fa- The mesonic loop correction part Vi 
becomes 

1 



Vi 



2K 



■[A^{T,fi) + 3A^{T,fi)] , 



(23) 



d'^p ( 1 



(24) 



where 



+ 



f d\ 1 
J (27r)3 2£;q 
d^q 1 



[1 - f{E^ + /.) - f{E^ - /.)] 



(27r)3 2E^+^ 



(25) 



Re/2(^a;^p) = -^ ^ 



Nf f d^q 1 



(27r)3 Ep+qEq 



^"^^ -.{f{Ep+, + ^^) + f{Ep+,-^^)-f{E, + ^^)-f{E,-^^)) 



" + El,, + ElY - AEl^^El 
x{f{E, + /.) - f{E, - /.) - f{E^, + ^,) + f{E^, - /.)) , (27) 

with o;^ = 27m/P and f{E) = l/{e^^ + 1)- the hmit of zero temperature, these equations 
are reduced to simpler forms: 

73, 



Vn 



T-»0 



^((M - + 0i) - 2^/ / [/^ + (^P - - /^)] ' (28) 



Ii{iu;,p) = 2Nf 



(29) 



(30) 



Rel2{iuj,p) = — / 



1 



X 



Nf 

4 J (27r)3Ep+qE, 

Ep+q — Eg 

a;2 + {Ep+g - EgY 



:{e{Eg-ii)-e{Ep+q-li)) 



Ep+q + Eq 



a;2 + (£;p+q + EqY 



{e{Eq-i,)-re{Ep+q-ij)) 



(31) 



Im/2(2CJ,p) 

'i (2^) 



a; 



{9{Eq - /.) - ^(£;p+^ - //)) , (32) 



where note that uj becomes a continuous valuable. 



In the present analysis, we adopt the parameter set of Hatsuda and Kunihiro [19| which 
are determined so as to reproduce F.^ = 93 MeV and = 138 MeV at T = /i = 0; the 
resultant values are rriQ = 5.5 MeV, the four-Fermi coupling constant G = 32.976 GeV^^ 
and the ultra-violet divergence cutoff yl = 631 MeV. 

Figure [1] represents the effective potential at /i = 330, 340, 348 and 500 MeV. When 
/i = 330 MeV, there appears a minimum around M = 340 MeV at the mean field level. 
The stability of the mean field solution can be investigated by the curvature of Vq in the 
(fa and ipT, directions. In the region denoted by "a unstable" , the curvature is negative and 
then unstable in the (p^ direction. Similarly, in the region denoted by "vr unstable", the 
curvature is negative in the (^,r direction. Fortunately, the minimum around M = 340 MeV 
is out of the unstable regions. This property is held for other values of fi; three examples are 
shown in other panels of Fig. [H As another interesting point, any vr unstable region does 
not appear for fi > 330 MeV. Thus, the mean field solution at the minimum point is stable 
at any fi for the case of the present parameter set. 

The meson loop corrections in the unstable regions do not make sense, since mesons 
considered there are tachyonic. Actually, the Gaussian integral in Eq.(j4]) breaks down for 
tachyonic mesons. Note that dashed curves in the unstable regions are just a guide of eyes. 
In the Nc = 3 case, we can not see where is a minimum, since it is somewhere in the unstable 
regions. So we take a somewhat larger N^, or 20, in which a minimum is still out of the 
unstable regions even after the inclusion of the correction. We then look into the effect of 
the next-to-leading order correction by comparing the mean field solution (the Nc ^ oo 
case) and the finite Nc (= 20) case. 

First we discuss the result of the mean field approximation. Studies in the mean field 

Q, [l4|. At /i = 340 MeV 
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approximation level have already been done in Refs. [11, 
(Fig. [11(b)) the minimum is still located around M = 340 MeV; this means that chiral 
symmetry is broken. At /i = 348 MeV (Fig. W^c)) two minima degenerate; in other words 
this is the first order phase transition point. At higher /i (Fig. [D^d)) chiral symmetry is 
restored to some extent. Here a comment is in order. In this case M at the minimum is still 
around 50 MeV, and it decreases gradually as fi increases; in other words we see a crossover 
to the restored phase. This situation was seen also in Refs. [ill . although not discussed. 



Now we consider the next-to-leading order correction due to a finite A^^^. In the case 
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Fig. 1: Non-dimensionalized effective potential as a function of the constituent quark mass. Cases 
of four different chemical potentials are presented. Thickness of curves decreases as the number 
of color increases as 3 ^ 20 ^ oo. Dashed curves indicate regions unstable with respect to the a 
and/or vr directions. For /x = 330 MeV, regions unstable in the a (vr) direction are denoted by "o" 
unstable" ("vr unstable"). For > 330 MeV, all regions denoted by dashed curves are unstable 
only in the a direction. 

of Nc = 20 the results are similar to the mean field case but the correction weakens the 
transition, in other words makes the jump in M small and shifts the critical /i to a lower 
value. These are clearly shown in Fig. [21 

To summarize, we have studied semi-quantitatively the effects of the mesonic, i.e., the 
next-to-leading order in the 1/Nc expansion, correction in the Nambu-Jona-Lasinio model 
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Fig. 2: The value of the constituent quark mass at which the effective potential becomes minimum 
as a function of the chemical potential. The mean field and the Nc = 20 cases are presented. 

on the high density chiral phase transition based on the auxiliary field method. The finite 
correction weakens the first order phase transition and shifts the critical chemical potential 
to a lower value. At Nc = 3, however, we can not see the minimum because of the instability 
of the mean field effective potential to the direction of the a and/or vr classical fields. For 
the mean field calculation, we have pointed out explicitly that chiral symmetry is not com- 
pletely restored immediately at the first order transition; after the first order transition the 
constituent quark mass decreases gradually — this may be called the two step restoration. 

We still have a lot to do: First of all we have to redetermine the parameter set, the 
four-Fermi coupling constant G and the ultra-violet cutoff A in the next-to-leading order. 
Next we have to study the phase diagram at finite T and fi to find the location of the critical 
end point. Further, studies of the three flavor model and the color superconductivity phase 
are to be done. 
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